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ON THE IMPEDANCE OF A SATELLITE BORNE 
VLF ELECTRIC FIELD ANTENNA 


Abstract: The source impedance of a satellite borne VLF antenna is calcu- 
lated using probe theory. Numerical examples of the impedance are given 
for typical plasma parameters for a short cylindrical antenna. It is shown 
that at these frequencies the resistive component dominates the imaginary 
component along both ionospheric and magnetospheric satellite orbits. This 
result is used to re-examine the interpretation in terms of electrostatic 
waves to explain the VLF noise observed by the satellite 1964-45A. It is 
shown that if a VLF antenna is monitored with a charge sensitive ampli- 
fier, the output signal is proportional to the product of the ambient signal 
level and the plasma density. The sustained noise enhancements observed 
on the 1964-45A satellite are thus attributed to changes in the antenna im- 
pedance along the satellite orbit rather than from changes in the ambient 


plasma noise level. 



Introduction 


The VLF electric field experiment aboard the 1964-45A satellite con- 
sisted of a short cylindrical antenna and four RMS voltmeters which mon- 
itored the AC voltages induced on the antenna in four frequency channels 
from 1.7-kc/s to 14.5-kc/s [Scarf et al., 1964]. The authors reported that 
background VLF electric field strength rarely fell below one mV/m in the 
ionosphere and the lower magnetosphere. In addition, sustained noise en- 
hancements were observed on the night side of the orbit which correlated 
with specific L shells and with the precipitation of energetic electrons. 

The experimenters interpreted these noise enhancements as direct evidence 
of electrostatic ion waves in those regions where the noise enhancements 
were observed. 

These experimental results are of considerable interest from the view 
point of plasma dynamics as well as geophysics. We have considered it 
worthwhile to consider in some detail the possibility that the observed 
noise enhancements were not field oscillations in the local plasma (i.e., 
geophysical phenomena) but might instead simply represent variations in 
the antenna impedance with changes in the ambient plasma parameters. 
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The Resistive Component of the Source Impedance 


Following the analysis of Mlodonsky and Garriott [l962] we will assume 
that at VLF frequencies the undisturbed plasma acts as an almost perfect 
conductor. In this case, which will be justified numerically later, the an- 
tenna impedance is determined primarily by the sheath impedance. Con- 
sider a short cylindrical antenna like that of the 1964-45A satellite as shown 
in Figure 1. The antenna will attain an equilibrium potential, 0, which is 
given by the solution of the integral equation 


I 


J e (0) • dS 



J + (<*>) • dS 



(4>) • dS + i 


(1) 


where J e (<£) is the electron current, J + (<£) is the ion current, J p (0) is the 
photo-emission current, and I is the current drawn from the antenna by the 
electrical loading of the receiver. Other charging effects such as secondary 
emission will be neglected. Equation (1) represents the condition of equilib- 
rium. This condition will be valid in time varying electric fields with fre- 
quencies less than the inverse of the relaxation constant r 1 of the plasma 
immersed antenna as will be discussed later. 

If a potential gradient exists over the dimensions of the antenna, or al- 
ternately if the antenna is moving in a magnetic field, the net current into a 
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small surface of the antenna will not necessarily be zero since the plasma 
potential will be a function of the position r of the surface element dS. In 
equilibrium, however, the total current to the antenna must again be zero. 



(r, 4 >) ■ dS 


!'• 


f 


(r, 0) • dS + T (r, 0) • dS + I 


( 2 ) 


If we assume that the electron distribution in the plasma is essentially 
Maxwellian and also that the potential of the antenna is negative, the elec- 
tron current can be treated analytically as 



Here 4> 0 is defined as the potential of the antenna with respect to the plasma 
potential at the origin of the coordinate system (see Figure la), J eo is the 
electron random current density, V is the velocity of the satellite, and B is 
the magnetic field strength. The electric field used to derive equation (3) 
was assumed to be constant over the dimensions of the antenna, i.e., the 
antenna length is taken to be short compared to the wavelength. 
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Combining equations (2) and (3) and integrating gives 


2tt J a KT d exp 


e(E + V x B) • d 


e(E + Vx B) • d 




J p • dS+I 


where a is the antenna radius and d is the antenna length. Equation (4) can 
be differentiated to yield the following expressions for the voltage gain and 
the source impedance of the antenna: 
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For short antennas e(E + V x B) • d < KT and the voltage gain can be approxi- 
mated by an expansion of the first term of equation (5) as 


E • d 
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The second term is a small correction due to focusing effects. When the potential 
is positive, the following simple empirical expression for the photoelectric current 
can be derived which agrees well with the experimental results of Hinteregger 
[l959] for energies below 1.5 volts: 


J P M 


Jpo eX P 


11-8 0 + 
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( 8 ) 


where J po is the photoelectric current density with no applied electric field 
and hv is a constant equal to 6.2 ev. 

Following the same procedure as in the case of negative potential, it is 
readily proved that 
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Equation (9) can be differentiated with respect to the electric field to 
yield the following expression for the voltage gain of the antenna: 

( 10 ) 


5< V E • d 

1 

hv 

d 

dE 

hv ^ 

jj e • dS-Jj + • dS 

<9E “ E 

I'll. 8 

exp j hv (E + Vx B) • dj 

|._ x 11.8(E+VxB) • d 

11.8 


*- d5 -. 

/• 

J+ ■ dS~I 


5 


11 8 -» -> -♦ -* 

Using the approximation (E + V x B) • d < 1, the above equation can be 
again approximated by an expansion of the first term on the right 


d4> 0 + 
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The source impedance is obtained by differentiating equation (9) with re- 
spect to current drawn from the antenna by the electrical loading of the 
receiver 



The voltage gain, equations (7) and (11) represents the response of an 
unloaded antenna with respect to the plasma potential at the origin to 
an external electric field. The actual AC voltage measured by a satellite- 
borne receiver will be larger than this amount because of the variation of 
the satellite ground potential with respect to the origin. This effect can be 
roughly approximated by replacing the antenna length d in equation (7) and (11) 
by an effective antenna length d e = d + t where t* represents the extent of 
the satellite along the direction of the vector d (see Figure la). Also it should 
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be pointed out that Equations (7) and (11) indicate that the voltage gain should 

— » — * 

be somewhat sensitive to the antenna orientation with respect to the vector V x B. 

Equations (6), (7), (11) and (12) are valid for time varying electric fields 
with periods longer than an equilibrium charging time r. This character- 
istic relaxation time is given roughly by the relation 

T = R S C S 

where C s is the antenna capacity which includes the sheath corrections. In 
addition equations (6), (7), (11), and (12) are strictly valid only at frequencies 
below the electron plasma frequency since the treatment given here ne- 
glects collective effects in the plasma and the effect introduced by the 
transit time of the electrons across the sheath. The treatment of such 
effects is beyond the scope of this discussion. This restriction does not 
limit the present investigation which is concerned with satellite borne VLF 
antennas since the plasma frequency in the magnetosphere is in general 
greater than 100 kc/s. 
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The Antenna Capacity 


It was suggested by Storey [l963] and Whale [l964] that the capacity of 
a plasma immersed cylindrical antenna might be approximated by that of 
the co-axial capacitor formed by the antenna and the boundary of the plasma 
sheath, i.e., 

C - 277 e 0 A. + a 

ln ^~~ 

where A. is the sheath thickness, a is the antenna radius, d is the antenna 
length, and e Q is the permittivity of free space. 

The above expression is not in general true since it diverges for 
zero potential. A more general expression can be derived by assuming 
a charge distribution in the sheath like that of Figure lb, where 
a + k + aA D is the distance at which the space charge becomes zero, A. D is 
the Debye length, and a. is a constant of the order of unity. 

The charge inside a cylinder of radius r and length d is given by 

q(r) = 277 d 



8 


where /j(r) is the charge distribution. The total charge can be calculated 
from the relation 
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a +A.+a\ r 


q(\) - 277 d I p(r) r dr 


(13) 


Using Gauss's theorem and E = -V^the potential on the antenna is obtained 


" 2tt 
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D [q(\)-q(r)] 


dr 


(14) 


The capacity of the antenna can conveniently be defined as 


dq(\) _ dq(MM 

d^o O) d0 o /dA. 


(15) 


Differentiating equation (14) and substituting into (15) we can obtain a simple 
expression for the capacity 
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In 
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(16) 


9 


Note that the above derivation does not require a knowledge of the exact 
nature of the charge distribution; only the distance k + aX. D is re- 
quired. The total charge inside the sheath can be found after integration of 
(13) between the limits a and a + k by assuming a constant charge density n e e 
shown in Figure lb. The result is 

q(k) = 7rdn e e [(\+a) 2 -a 2 ] (17) 

Combining equations (14) and (17) we can find the value of k from the 
expression 

n e e f a + \ + ak D 1 'I 

4> 0 = 2e “ |( a + M 2 In 2 [( a + X + a \>) 2 “ a 2 ]| ( 18 ) 

when the value of 4> 0 is known. 

The above treatment of antenna capacity and sheath thickness does not 
include the effects of the satellite wake and of the antenna ground plane. 

It should be noted that since the plasma parameters enter only into the 
log term of the antenna capacity the antenna capacity should not vary ap- 
preciably as a function of the satellite plasma environment. In 
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contra-distinction the resistive impedance of the antenna is a much stronger 




function of the plasma parameters varying roughly as R s ~ 


KT 

n e 


Numerical Calculations 

The ion current to the antenna when the potential on it is negative can 
be calculated using the following expression, [orsak et. al., 1965] 

/ KT \ 1//2 

J+ M = e \2Fsrj n + Wcr (i9) 

where: ¥ is an ion current factor due to the accelerating potential 4> 0 (see 
appendix) and cr is an ion current factor due to satellite velocity (see appendix). 

To calculate the factor, the floating potential is needed. An approx- 
imate value of it is given by 




KT 1 eo 

e 1° T 
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where I is the electron random current and I. _ is the ion ram current. 

eo “ k 


li 


The sheath thickness can be obtained from equation (18), using this cal- 


culated value of 4> 0 . Inserting <t> 0 and A. in the equation (A-l), the V factor is 
determined. 

Thecr factor can be calculated from the known ion velocity, if the sat- 
ellite velocity is given. For roughly circular orbits 

0.79 • 10 6 

v - yiTR^ cm/sec 


where 

R = radius vector 
R 0 = earth radius 

Wave effects will reduce the ion current somewhat. 

The photo emission current can be estimated from the experimental 
results of Hinteregger [l959] : 
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po 


• dS 


3.9x 10 9 S is [Negative potential] 


( 20 ) 
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where 


S is = 2ad for the cylindrical antenna. 


Shadowing will of course reduce this current. Since the photoelectric current 
does not change appreciably with altitude as long as the potential on the an- 
tenna remains negative and since the electron and ion current decreases 
with increasing altitude, it is expected that the potential on the floating 
antenna will become zero in a specific altitude. This occurs at R/R 0 = 2.6 
which can be derived by setting <£ 0 = 0 in equation (1) and using equations 
(19) and (20) for the ion and photoelectric current respectively. 

In the case of positive potential the electron current is given by 

/ KT \ 1/2 

J.(v) = '(2^ r) ■>.* 

The factor, similar to the previous case, is due to the accelerating po- 
tential (see Appendix). In this case, however, the contribution from the 
motion of satellite can be neglected since the electron thermal velocity is 
much greater than the satellite velocity. 
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While the ion current could be calculated from the expression [Kanal, 
1964] 



(here y 2 - 


V 2 m. 

2KT 


, V is the satellite velocity and I is the modified Bessel 


function of the n th order), this is not necessary for the most part 
since its maximum value is only a few percent of the photoelectron 
current and therefore can be omitted without introducing a considerable 


error. 


The value of <£ 0 + which is necessary for the computation of the *P factor 
can be approximated by the expression 
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while the sheath thickness is given again from equation (18) 
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In order to give numerical examples of the antenna impedance, typical 
values of the plasma ambient plasma parameters n and T are needed. A 
useful graphical summary of the plasma conditions in the ionosphere as 
given by Bourdeau [l965] and Johnson [l962] are shown in Figures 2a and 2b. Es- 
timates of the thermal plasma parameters in the magnetosphere are slightly 
less well founded at present. Figure 3 shows an average of several orbits 
of the ion trap measurements of IMP-II [Serbu and Maier, 1966] and one 
orbit of IMP-I measurements. Calculated values of the resistive portion of 
the source impedance of an antenna of 1 cm diameter and 100 cm length as 
a function of altitude are shown in Figure 4 for the ionosphere and Figure 5 
for the magnetosphere. These graphs were calculated using equation (6) 
for negative potential, equation (12) for positive potentials, and the plasma 
parameters given in Figures (2) and (3). The axis of the antenna was taken 
as perpendicular to the solar vector, perpendicular to the satellite velocity 
vector, and perpendicular to the V x B vector. It can be seen from these graphs 
that the resistive portion of the impedance varies considerably with alti- 
tude. Values for the geometry chosen vary from Kilo-ohms in the F-layer 
to tens of Meg-ohms in the outer magnetosphere. 

Values of the capacity of an antenna of similar dimensions were calcu- 
lated from equation (16) and are shown as a function of altitude for the 
ionosphere in Figure (6) and for the magnetosphere in Figure (7). The 
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antenna capacity is much less dependent on altitude varying from about 9 pf 
minimum in the outer magnetosphere to 55 pf maximum in the daytime 
F -layer. 

The above calculations for the antenna impedance are appropriate to a 
satellite-borne antenna. The resistive component of an antenna aboard a 
sounding rocket near apogee will be an order of magnitude larger since the 
ion ram current will be much less. It should also be pointed out that these 
calculations are for an isolated antenna at the floating potential. By biasing 
electron current from the antenna into the satellite the resistive impedance 
can be considerably reduced [storey, 1963]. 

It is important to compare the resistive source impedance with the 
capacity impedance. Since these impedances act electrically in parallel on 
the antenna [Mlodonsky and Garriott, 1962] the resistive impedance will 
dominate at frequencies below a value f 0 given by the expression 


R 1 

“• ~ 2nf 0 C s 

Values of this transition frequency f 0 are given in Figure 8 for the 
ionosphere and in Figure 9 for the magnetosphere. It can be seen from 
these curves that the resistive impedance of a satellite-borne antenna 
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dominates the capacitive impedance at VLF frequencies in both the iono- 
sphere and in the lower portions magnetosphere. 

The transition frequency f Q is also equal to the upper limit of the 
validity of the derivation of the resistive impedance, equation (6) and (12), 
since above this frequency the antenna cannot reach the assumed equilibrium 
condition. 

It is important to justify the neglect of the bulk impedance of the plasma 
if we wish to apply the calculated sheath impedances directly to the in- 
terpretation of experimental data. The bulk conductivity of the plasma is 
extremely anisotropic due to the magnetic field being a fair conductor along 
the field vector and a good insulator perpendicular to it. The conductivity 
along the field vector is given by the Lorentz conductivity 

V - ICO 

m 


where is the plasma frequency and v m is the collision frequency for 
momentum transfer of the electrons. Since the collision frequency is con- 
siderably less than the signal frequency at VLF frequencies along most 
satellite orbits, the bulk conductivity as given above is mostly inductive 
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[Storey, 1963], One can estimate the frequency domain where the sheath im- 
pedance dominates the bulk impedance as follows 



co C 
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<< 


oj 2 e 0 d 


CO 


(21) 


For VLF frequencies it has been numerically shown that 


IT >> -C s (22) 

s 


Using this result and the approximation C s - e 0 d equation (21) reduces to 

i rk~ <K < (23) 

S S K 

From (22) and (23) we have a frequency domain where the sheath impedance 
dominates the bulk impedance 


CO « CO 
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(24) 


18 



This frequency domain includes the VLF phenomena at satellite altitudes as 
already noted. 

Antenna Pre-Amplifier Considerations 

From the above considerations we can conclude that although the volt- 
age gain, equations (7) and (11), of a satellite-borne VLF antenna is to a 
first approximation independent of the ambient plasma parameters, the an- 
tenna impedance is a strong function of these parameters. Of the two 
satellite VLF electric field intensity experiments reported to date. The 
experiment of Storey [l963] utilizes a voltage sensitive preamplifier as 
shown schematically in Figure 10a. Modern semi-conductor technology 
allows the design of such an amplifier with a resistive input impedance of 
the order of hundreds of Meg-ohms. Thus this type of preamplifier allows 
the measurement of VLF signals with receiver gain independent of the 
ambient plasma conditions. 

On the other hand, the VLF experiment on the 1964-45A satellite utilized 
a charge sensitive pre-amplifier [scarf et. al., 1965] as is illustrated 
schematically in Figure 10b neglecting the sheath capacity. The output 
signal of such an amplifier is proportional to the time integral of the input 
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(25) 
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ou t 



dt 


where V s would be the signal induced on the antenna if it was unloaded by the 
amplifier and R s is the sheath impedance. At a particular frequency f n in the 
Fourier domain this expression reduces to the form 



(V.) 


fn 


R 


at the altitudes of this experiment (300 km - 4000 km) the source impedance 
is determined primarily by the ion ram current 


KT 1 

R s a ~ en VS,' 

Thus, to a first approximation, this experiment apparently measured a 
quantity proportional to the product of the ambient electric field intensity 
and the plasma density 

E f n e 

V out a KT ( 26 ) 
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Alternate Interpretation of the Measurements of Scarf et al. 


In view of the above analysis of the actual physical quantity measured in 
the 1964-45A satellite electric field experiment we feel that a reinterpre- 
tation of the geophysics involved is in order. As mentioned earlier the RMS 
voltages measured indicated a slowly varying background level on the day- 
time apogee half of the orbits with sustained signal enhancements at 1.7-kc/s 
on the nighttime perigee portion with some general correlation with specific 
L shells and with the precipitation of energetic electrons. The authors 
[Scarf et. al., 1965] interpreted these signal enhancements as direct evidence 
of electrostatic ion waves in the ambient plasma since they could not conceive 
of an alternate explanation of the correlation of the signal enhancement with 
specific L shells and energetic electron fluxes. 

The point to be made here is that there is an alternate explanation of 
the observed correlation of these signal enhancements with either specific 
L shells or energetic particles in terms of the variation of the antenna- 
amplifier response with the ambient plasma density as given by equation (26). 
The magnetic control of the plasma density in the upper ionosphere was ob- 
served with a probe aboard the Ariel I satellite [Sayers et. al., 1963] and 
confirmed with the first top side sounder satellite [King et. al., 1964]. In 
addition, observations with both the Ariel and Aloutte I satellite suggest that 
energetic particles should be considered as an ionization source of the 
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F-region electron density [Bordeau, 1965]. More recently Sharp [l966] has 
found troughs and enhancements in the plasma density as a function of latitude 
by means of a ion trap carried aboard an earth-oriented circular satellite 
at roughly 300 km altitude. The total variation in plasma density during one 
orbit was larger than a factor of one hundred. Converting these changes in 
plasma density to changes in antenna amplifier response via equation (26) 
allows one a simple explanation of the apparent noise enhancement measured 
aboard the 1964-45A satellite in terms of changes in antenna impedance. 

In addition to the variations in the relative noise levels in the ionospheric 
plasma, the absolute magnitude of the noise is of great interest. Since the 
analysis given in this paper indicates that the 1964-45A electric field ex- 
periment monitored a quantity which was physically different from what it 
was calibrated for, we feel that the absolute magnitudes of the noise levels 
reported are perhaps also questionable. 

Conclusions 

I. The resistive component of the source impedance of a satellite-borne 
antenna dominates the capacitive component at VLF frequencies in both the 
ionosphere and the lower magnetosphere. 

II. The antenna impedance varies considerably as a function of altitude. 
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III. The experimental results of the VLF electric field experiment [Scarf 
et. al., 1965] flown aboard the 1964-45A satellite can perhaps be re- interpreted. 
In particular the correlation of the apparent noise enhancements observed 
with specific L shells and with the observation of precipitating energetic 
electrons can be interpreted as representing changes in the antenna im- 
pedance because of the charge sensitive characteristics of the amplifier 
used rather than representing an observation of ion-acoustic waves. 
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appendix 

The ion current density to the antenna, for negative potential and assum- 
ing thermal equilibrium of ions, is given by [orsak et. al., 1965] 


where : 
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ion random current 


n + = ion density 
m. = ion mass 
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a = antenna radius 
A. = sheath thickness 
T = ion temperature 
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The a factor represents the contribution to ion current from the satellite 
motion. Assuming that the velocity vector is perpendicular to the antenna 
axis it can be expressed as 


Cr 



where 




4>{ x ) dx - 1 


0(x) 



and 


x 



1/2 


and E is a complete elliptic integral of the second kind. 
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FIGURE CAPTIONS 


Figure la 
Figure lb 
Figure 2a 


Figure 2b 


Figure 3 


Figure 4 


Figure 5 


Antenna geometry 

Postulated charge distribution inside the sheath 

Typical electron density for the ionosphere as summarized 
by Bordeau [l965] 

Typical electron temperature for the ionosphere [Bordeau, 
1965] 

Plasma parameters for the magnetosphere [Serbu and 
Maier, 1966] 

Resistive impedance of a short cylindrical antenna (0.5 
radius; 100 cm length) in the ionosphere. The antenna axis 
assumed perpendicular to the direction of motion and 
T = T. , n = n. . 

Resistive impedance of a short cylindrical antenna (0.5 
radius; 100 cm length) in magnetosphere. The antenna axis 
assumed perpendicular to velocity vector and the solar 
vector. Also T = T. and n = n. . 
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Figure 6 


Figure 7 


Figure 8 


Figure 9 


Figure 10 


Capacity of a short cylindrical antenna (0.5 radius; 100 
cm length) in the ionosphere. The antenna axis assumed 
perpendicular to the direction of motion. The cut-off 
parameter a is taken as unity for these calculations. 

Capacity of a short cylindrical antenna (0.5 cm radius; 100 
cm length) in magnetosphere. The velocity vector assumed 
perpendicular to the antenna axis and the solar vector and 
the cut-off parameter a is taken as unity. 

Transition frequency for a satellite-borne VLF antenna in 
the ionosphere. The region to the left of curves represents 
the domain where the resistive part of the impedance domi- 
nates the capacitive part of the impedance. 

Transition frequency for a satellite-borne VLF antenna in 
the magnetosphere. The region to the left of the curve 
represents the domain where R. dominates of R c . 

Preamplifier schematics. V. represents the unloaded source 
voltage and R s represents the sheath impedance. The ca- 
pacitor in the charge sensitive amplifier may be alternately 
returned to ground rather than the output terminal with no 
appreciable change in the transfer functions over the fre- 
quency region where the amplifier is charge sensitive, R s C> 
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Figure 7 
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a) Voltage Sensitive Amplifier 
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